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Summary 

An analysis has been performed to determine the 
effects of model error on the control of a large flexi- 
ble space antenna. Control was achieved by employ- 
ing two three-axis control-moment gyros (CMG’s) lo- 
cated on the antenna column. State variables were 
estimated by including an observer in the control 
loop that used attitude and attitude-rate sensors on 
the column. Errors were assumed to exist in the in- 
dividual model parameters: modal frequency, modal 
damping, mode slope, and moment of inertia. Their 
effects on control-system performance were analyzed 
either for (1) nulling initial disturbances in the rigid- 
body modes, or (2) nulling initial disturbances in 
the first three flexible modes. The study includes 
the effects on stability, time to null, and control re- 
quirements (defined as maximum torque and total 
momentum), as well as on the accuracy of obtain- 
ing initial estimates of the disturbances. The effects 
on the transients of the undisturbed modes are also 
included. The results, which are compared for de- 
coupled and linear quadratic regulator (LQR) con- 
trol procedures, are shown in tabular form, para- 
metric plots, and as sample time histories of modal- 
amplitude and control responses. 

Results of the analysis showed that the effects of 
model errors on the control-system performance were 
generally comparable for both control procedures. 
The effect of mode-slope error was the most serious 
of all model errors. 

Introduction 

Large space structures (LSS) will require control 
for attitude orientation as well as procedures for con- 
trolling the flexible modes. After deployment or con- 
struction, any control required will most likely be 
performed with a control procedure based on inac- 
curate knowledge of model parameters. There are 
many identification and adaptive control procedures 
that have been developed (for example, ref. 1) for 
updating and verifying these parameters. However, 
the practical application of these procedures for op- 
eration in space has not been fully tested. Hence, a 
considerable period of time may be required before 
a more accurate model can be determined. Conse- 
quently, a knowledge of the effects of model errors on 
the stability and control requirements for LSS control 
would be useful. 

The effects of model errors are addressed in 
the present paper with regard to the control of a 
122-m-diameter hoop-column antenna (ref. 2). A 
finite-element mathematical model of the antenna 
is used that provides a convenient tool for simu- 
lated control analysis of realistic large flexible space 
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Solar panels 


Sketch A 

structures. Previous control studies that used this 
antenna model include references 3 to 5. The an- 
tenna is depicted in sketch A and is described in 
reference 2. 

The present paper is based on the control analysis 
reported in reference 4. In this reference an observer 
was used for estimating the state variables for control 
feedback. In addition, residual modes were included 
in the antenna model and the effects of including 
these modes were investigated. The control arrange- 
ment employed one three-axis control-moment gyro 
(CMG) at the top of the antenna column and one at 
the bottom. Only the first three flexible modes were 
included in the control law, inasmuch as data in ref- 
erence 3 showed large increases in the control require- 
ments with additional modes included. In reference 4 
it was shown that decoupled and linear quadratic reg- 
ulator (LQR) control procedures gave comparable re- 
sults; hence, the effects of model errors for both pro- 
cedures are analyzed and compared in the present 
analysis. 

Errors are assumed to exist in the individual 
model parameters: modal frequency, modal damp- 
ing, mode slope (control-influence coefficient), and 
moment of inertia. Their effects on control-system 
performance are analyzed either for (1) nulling initial 
disturbances in the rigid-body modes, or (2) nulling 
initial disturbances in the first three flexible modes. 
The study includes the effects on stability, time to 
null, and control requirements (defined as maximum 
torque and total momentum), as well as on the ac- 
curacy of obtaining initial estimates of the distur- 
bances. The effects on the transients of the undis- 
turbed modes are also included. Results are shown 
in tabular form, parametric plots, and as sample time 
histories of modal-amplitude and control responses. 
For the analysis, random positive and negative errors 



were incorporated into the frequency and mode-slope 
parameters, as were certain errors in the damping 
and moment of inertia; their effects were analyzed 
individually for specified magnitudes of error. The 
analyses were carried out for a large number of sam- 
ples in order to obtain the maximum effects of these 
errors. The control system and observer were as- 
sumed to be perfect in that no actuator or sensor 
dynamics were included in the analysis. 


Symbols and Abbreviations 


A 

system matrix (eq. (4)) 

An 

modal amplitude (eq. (2)), 
where n is mode number 
(1,2,... ,6) 

B 

control influence matrix 
(eq. (4)) 

C 

observation matrix (eq. (9)) 

CMG 

control-moment gyro 

d 

displacement vector (eq. (3)) 

F 

decoupled feedback gain 
matrix (eq. (5)) 

G 

decoupled feedforward gain 
matrix (eq. (5)) 

I 

moment-of-inertia matrix 
(eq. (1)) 

/ 

model moment of inertia 

Ixi lyi Izi Ixz 

center-of-gravity moments and 
product of inertia 

J 

objective function (eq. (7)) 

K 

LQR feedback gain matrix 
(eq. (6)) 

K' 

estimator gain matrix 
(eqs. (13)) 

LQR 

linear quadratic regulator 

l 

distance (see fig. 1) 

M% ,j 

angular momentum, where i 
is control direction and j is 
actuator location 

m n 

nominal total momentum (no 
model error) 

M e 

total momentum with model 
error 


Q 

state-vector weighting matrix 
(eq. (7)) 

R 

control-vector weighting 
matrix (eq. (7)) 

r 

radius of hoop 

s 

vector of system input noise 
(eq. (10)) 

T 

torque vector 

T ■ 
1 i,3 

components of torque, where 
i is control direction and j is 
actuator location 

Tmax 

maximum value for actuator 
requiring largest torque 

T N 

nominal |T max | (no model 
error) 

T n 

natural period 

T e 

|T max | with model error 

t 

time 

U 

vector of observation noise 
(eq. (8)) 

V 

input command vector 
(eq. (5)) 

w 

weight 

X 

state vector 

X 

estimate of controlled state 
vector 

x, y, z 

coordinates of antenna center 
of gravity 

y 

observation vector (eq. (8)) 

z 

sensor output vector without 
noise 

e 

estimator error, y — C'x 

& 

desired damping ratio 

$n 

natural damping ratio 


rotation angle about a;-, y-, 
and 2-axis, respectively 

$ 

mode-shape matrix (eq. (3)) 


mode-slope matrix, where n 
is mode number, i is control 
direction, and j is actuator 
location 


desired frequency 


natural frequency 
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Subscripts: 

A n modal amplitude for nth mode 

i , j control direction and actuator 

location, respectively 

0, 0, $ rotation angle about x-, y- y 

and 2-axis, respectively 

A matrix with a prime relates to the estimator 
equations. Dots over symbols indicate derivatives 
with respect to time. A circumflex 0 over a symbol 
indicates an estimate of the state variable. All 
matrices are given in units of inches, pounds, and 
seconds. 

Mathematical Model of Antenna 

A large space structure such as the hoop-column 
antenna has, in theory, an infinite number of flexible 
(vibration) modes. To facilitate analytical treatment 
of the control problem, a finite-order linearized model 
was formulated. For the analysis of this report, the 
structural model was selected to contain six of the 
lowest flexible modes of the 122-m-diameter hoop- 
column antenna, as described in reference 2. The 
three rigid-body rotation modes are included in the 
analysis. 

The equations of motion used to represent the 
rigid-body and vibration modes of the antenna are 
now discussed. Rigid-body rotations (for small an- 
gles) about the antenna center of gravity for two 
three-axis CMG actuators on the column are rep- 
resented by 



where 


T = [Tx,u Ty, 1, T zA , T x>2 , Ty 2, T z $ 


equation (2) are modal-amplitude displacement vari- 
ables and do not represent actual physical displace- 
ments. The physical displacement at some point on 
the antenna is a linear combination of the modal 
amplitudes and mode shapes and is given by the 
transformation 

d j = <&A n (3) 

where <& is the mode-shape matrix. 

The mode-slope data used in the analysis were 
taken from results of a NASTRAN® model of the 
122-m-diameter hoop-column antenna system given 
in reference 3. Table I shows the weight and inertias 
of the antenna. The frequencies and assumed natural 
damping ratios of the six flexible modes considered 
in the analysis are also given in table I. As a matter 
of interest, the unforced transient behavior of these 
modes is presented in figure 2. 

Decoupled and LQR Control 

The second-order equations in the analysis 
(eqs. (1) and (2)) can be reduced to first-order equa- 
tions (state-vector form) and written as 

x = Ax + BT (4) 

where the state vector 

x — [0, <t>, Ai, A 2 , A 6 , 0, (j), Ai, A 2 , Aq} t 

Ai, A 2 , and A3 are the controlled flexible modes, and 
A4, A5, and Aq are considered as residual modes. 
The B matrix used in the present analysis is given in 
table AI in the appendix. The composition of the A 
matrix is defined as follows: 

°9x9 | 1-9x9 

Ai8xl8= - 

“^9x9 | “D9 x9 


The CMG’s were used on the column such that 
control torques were about the x— ,y— , and 2— axes 
as defined in figure 1. The CMG’s are numbered to 
be consistent with reference 3. 

Variations in modal amplitudes of the flexible 
modes are represented by 

An "T 2 $n^nAn T OJ n A n = 

m n 

(n= 1,2, ..,6) (2) 

where m n is the modal mass and is the mode- 
slope matrix. It should be noted that A n values in 


where the subscripts indicate the dimensions of the 
various elements, and 

0 null matrix 

1 identity matrix 

A diag(0, 0, 0, coj, a|) 

D diag(0, 0, 0, 2^2, 2&we) 

Decoupled Control 

Since decoupling theory was analyzed extensively 
in reference 3, the discussion of its application in 
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this paper is limited to some general remarks. The 
decoupling control law is taken as 


T = Fx + Gv (5) 


where x is the estimated state vector, v is the input 
command vector, and F and G are feedback and feed- 
forward gain matrices, respectively. (As noted, esti- 
mates of the modal variables are employed because 
these variables cannot be measured directly.) De- 
coupling theory provides a method for determining 
the F and G matrices such that independent control 
is maintained for each of the decoupled (controlled) 
state variables. Different values can be selected for 
the closed-loop dynamics (t and q) without affect- 
ing the independent control capability. Examples 
are given in reference 3 that illustrate the manner 
in which the theory was applied to the system equa- 
tions in determining the F and G matrices. The 
values for the decoupled closed-loop dynamics used 
in the current analysis are given in table II. This case 
corresponds to Q values that give approximately the 
same value (10 sec) for time to damp to 1 percent for 
each of the three flexible modes. The other case rep- 
resents the closed-loop dynamics pertaining to the 
LQR control analysis, which are close to those for 
decoupled control. The corresponding F gain matrix 
used in the current analysis for decoupled control is 
given in table All in the appendix. The G matrix, 
which is required for rigid-body commands, was not 
used in the present analysis. Rigid-body commands 
were not analyzed inasmuch as the results would be 
similar to those for the control for the rigid-body 
disturbances. The F matrix was calculated for a 
reduced-order model (no residual modes). Hence, 


x = 


0, 0, Au M, M, 0, 4>, Vh My 


M, 



LQR Control 

The asymptotic linear quadratic regulator (LQR) 
approach of reference 6 was also used to calculate 
the state feedback control gains. (A computer pro- 
gram is given in ref. 7 for performing these compu- 
tations.) When comparing the LQR and decoupled- 
control procedures, it should be noted that the de- 
coupled approach has advantages in some areas. For 
example, the decoupled feedback gains are computed 
in closed form and yield the exact desired closed-loop 
dynamics, whereas the LQR approach requires an 
inexact, time-consuming iterative solution to obtain 
approximate desired closed-loop dynamics. Also, the 
decoupled gains give a known transient response (for 
perfect state- variable estimates), as well as a steady- 
state response. Only the steady-state response is 
known by use of the LQR method. 

The control law for the LQR results was given by 


T — Kx (6) 

The control gain matrix K was computed subject to 
the constraint of equation (4) such that the following 
function was minimized: 

J = ^lim> jf* (x t Qx + T t RT) dt (7) 

where Q and R are positive definite symmetric 
weighting matrices for the states and controls, which 
can be varied in order to achieve desired closed-loop 
dynamics. For the present analysis, the identity ma- 
trix used for R was found to be sufficient; the state- 
variable weighting matrix Q was adjusted to obtain a 
set of closed-loop dynamics close to that used for de- 
coupled control. (See table II.) The corresponding K 
gains were used in the control procedure (eq. (6)) for 
the subsequent LQR control analysis. The weight- 
ings Q and the K-gain matrix are given in table AIII 
in the appendix. 


Incorporation of Observer 

An observer was incorporated into the control procedures in order to calculate estimates of the state variables 
as required in equations (5) and (6). The Kalman-Bucy filter, which represents a systematic procedure for 
observer design, is used for obtaining these estimates. An attitude-measurement device (e.g., a star tracker) 
and a rate-measurement device (e.g., a rate gyro) were assumed to be present at certain locations on the 
column. The sensor output is given by 


y = z + u 


( 8 ) 
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where z = Cx and u is the sensor noise, which is assumed to be 0 for this analysis. As an example, the 
observation matrix C is 


^6x18 = 






|* 2 ? * 5 ? • 

. *'(«) 
X,1 



03x9 


Ux3|<? <?■ 






* 2 ? * 2 .?- 

■* 2 6 ? 

1 

1 

# ' ( i ) 

x,4 

*' ( 2 ) . 

x,4 

■•*22 | 



1-3x3 | 

I 

*' ( i ) 

yA 

*' <2 > • 
yA 

■•«£? i 

03x9 


1 

1 

z,4 

■ 

z,4 

■•*22 i 




( 9 ) 


and 


x — [o, <f> t V, ^i, A 2 , A 6 , 9, j>, Ai, A 2 , .... A 6 ] 


The matrix I 3 X 3 * s the identity matrix and 03 x 9 is the null matrix. The superscripts on the mode slopes refer 
to the flexible-mode number, and the subscripts refer to the direction and location. For example, the subscript 
x, 1 refers to the mode slope in the ^-direction at location 1 on the antenna column (fig. 1.) The observation 
matrix shown in equation (9) corresponds to a sensed attitude rate in the x-,y-, and ^-directions at location 1 
and to a sensed attitude in the x - , y-, and ^-directions at location 4. Mode-slope values for various locations 
on the column are given in reference 3. The values corresponding to equation (9) are given in table AIV in the 
appendix. 

The complete set of equations required for the control analysis becomes 


x = Ax + BT + s 

where s, the process noise, is assumed to be 0 for this analysis, 


( 10 ) 


T = Kx (for LQR control) 

T = Fx + Gv (for decoupled control) 
y = Cx 


(lla) 

(lib) 

( 12 ) 


* = A'x + B'T + K'(y - C'x) 


where 


£ = (A' - K'C' + B'K)x + K'Cx (for LQR control) (13a) 

x — (A' — + B r F)x + I^Cx + B^v (for decoupled control) (13b) 

The estimator equations (eqs. (13)) calculate the estimates of the state vector using the measurement y and 
the estimator gains K 7 , which are determined by a steady-state Kalman-Bucy filter. If all six flexible modes 
are estimated, then C' = C. However, if the residual modes are not estimated, as in the present analysis, C f 
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becomes a 6 X 12 matrix. In this case, the six columns that contain the information are 

eliminated. The estimator gains K 7 were calculated by the technique given in reference 7. In conjunction with 
the C' matrix, it should be emphasized that there are two differences between the A', B', and C' matrices 
and the A, B, and C matrices. First, the primed matrices are of lower order than the unprimed matrices 
since they do not include the residual modes. Second, the unprimed matrices represent the actual (true) model 
parameters as revised by the various model parameter errors that modified the nominally assumed values. 

The set of composite equations used for the LQR-control analysis is 


/ X \ 

A 18 XI 8 ! 

1 

B 18x6K6xl2 

V x y 

K'i2x6C6xl8 I 
L i 

A'i2xl2 - K'^xeC'ex^ + B, 12x6K6xl2 


For the decoupled-control analysis, K is replaced by F in equation (14) and the term 



(14) 


B 18x6 G 6 x 6 
B'i2x6 G6 x6 


is added to the right side of the equation. The subscripts indicate the dimensions of the various elements for the 
case in which the three residual modes are not estimated, which was used throughout the current analysis. With 
precalculated values for K' and K (or K', F, and G), the composite set of equations was integrated to obtain 
time-history responses for specified initial conditions in the various modes. In determining the estimator gain 
matrix, the weighting matrices S and U for the states and observations, which are analogous to the covariance 
intensities of s and u, were varied until acceptable performance (satisfactory estimator eigenvalues and realistic 
gain values) was obtained. The resulting matrices were 

S = diag(.l, .1, .1, 10 5 , 10 5 , 10 5 , .1, .1, .1) 

U = diag( 1, 1, 1, 1, 1, 1000) 

The weighting value of 1000 gave the best results as far as control oscillations were concerned. 

A block diagram for decoupled control is shown in sketch B. (A block diagram for LQR control is the same, 
with the exclusion of the G matrix.) 



^Decoupled controller 


Sketch B 
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Method of Analysis 

For the study, errors were incorporated in the individual sets of model parameters, one parameter at a 
time. For example, the effect of errors in knowledge of the flexible-mode natural frequencies oj n was studied 
with the assumption that the values of all other model parameters were correct. Random sets of positive and 
negative errors were employed for investigating the effects of frequency and mode-slope error. (For example, see 
table III.) Effects of errors in damping and moment of inertia were investigated for either positive or negative 
errors. As previously described, the frequency errors are reflected in the plant matrix A, which is a top element 
in equation (14). Errors in the damping ratio are also associated with the A matrix. Errors in the moments 
of inertia are associated with the B matrix, and for errors in the mode slopes, both the B and C matrices 
are affected. The feedback gain matrices K or F, as well as the estimator gain matrix K 7 , were computed 
using the nominal (no errors) A', B', and C' matrices. Also, these matrices were employed for the estimator 
in equation (14). The resulting eigenvalues for the composite set of equations (see eq. (14)) were checked for 
stability for each set of random errors used. 


Results and Discussion 

Results of the study are presented to illustrate 
effects of various model errors on the control sys- 
tem performance with regard to stability, time to 
null, and control requirements. The effects are an- 
alyzed either for (1) nulling initial disturbances in 
the rigid-body modes, or (2) nulling initial distur- 
bances in the first three flexible modes. Case (1) 
is referred to as “rigid-body control,” even though 
control is being applied to attempt to maintain the 
first three flexible modes at 0. Similarly, case (2) 
is referred to as “flexible-mode control.” Compar- 
isons are made between the decoupled and LQR con- 
trol procedures with the results shown in tabular 
form, parametric plots, and as sample time histories 
of modal-amplitude and control responses. Compar- 
isons are also made for the effects on the transients 
of the undisturbed modes. 

Momentum and the maximum value for torque 
input are used as the measure of control requirements 
when making various comparisons. The momentum 
time histories represent the area under the torque 
response curves. The momentum values in the tables 
and parametric plots are the values at the end of 
each computer run, when the torque responses have 
essentially zeroed out. Maximum torque is the value 
for that actuator requiring the highest absolute value 
of torque. 

General Considerations 

For simulation, two three-axis CMG’s were em- 
ployed as actuators to control initial disturbances of 
0.01 rad in the rigid-body modes or 1 in. in the three 
flexible mode amplitudes. Since the system equations 
are linear, the results can be scaled to any constant 
positive multiple of these initial conditions. Initial 
disturbances were also included in the three residual 
modes, even though it was determined that the re- 
sults were not materially affected by their inclusion. 


(See ref. 4.) The CMG’s were placed at the top and 
bottom of the antenna column at locations 1 and 2 
as shown in figure 1. It was assumed that three- axis 
attitude sensors and three-axis attitude-rate sensors 
were available at various locations on the column for 
determining estimates of the state variables. Unless 
otherwise specified, one three-axis attitude-rate sen- 
sor at location 1 and one three-axis attitude sensor 
at location 4 were used. 

It was further assumed that the control system 
was not activated until the initial state estimates 
were known to within 10 percent of true values for 
the rigid-body modes and for the first three flexible 
modes. The initial state estimates are actually the 
initial conditions in the estimator equations. Esti- 
mates of the residual modes were not included in the 
control procedure. Further discussion concerning ini- 
tial estimates is given in the following section. 

Effect of Model Error on Initial Estimates 

As shown in reference 4, good estimates of the 
state variables are required before activating the con- 
trol system, especially for rigid-body control. With- 
out good estimates, the control requirements are ex- 
cessive. Good estimates can be assured by running 
the estimator a few seconds before turning on the 
controller. The time histories of estimator error for 
rigid-body control are shown in figure 3. As ex- 
pected, these time histories are unaffected by any 
type of model error up to the magnitudes considered 
in this study. The curves in figure 3 show that good 
estimates (90 percent) of the rigid-body disturbances 
are obtained within about 12 sec. 

Figure 4 shows examples of the important effects 
of model errors in determining the state estimates for 
flexible-mode control before the controller is turned 
on. The effects can be compared to the case with no 
model errors, as shown in figure 4(a). Figure 4(b) 
shows that with <&' errors of ±20 percent, state 
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estimation accuracy after about 6 sec is about 80 per- 
cent. It was determined that for errors of ±50 per- 
cent, the state estimates attained only about 50 per- 
cent accuracy in 6 sec. Figures 4(c) and 4(d) indicate 
that oj n errors are more crucial than errors with 
regard to the initial accuracy of the state estimates. 
For oj u errors of ±20 percent, the initial state esti- 
mates are shown to be good to only about 50 percent. 
For large ou n errors (±50 percent), the observer was 
marginally stable. Errors in the moments of inertia 
and damping ratios had practically no effect on the 
initial state estimates. 

As previously stated, the results in this paper per- 
tain to the assumption that the initial state estimates 
are within 10 percent of the true values at the time 
that the control is turned on. These results, how- 
ever, are applicable to larger percentage errors in the 
estimates, inasmuch as reference 4 shows that satis- 
factory control performance for flexible-mode control 
is obtained for initial state estimates that were within 
40 percent of the true values. 

Effect of Model Frequency Error on Control 

Performance 

Flexible-mode control . Figure 5 shows sample 
time histories of modal responses and control require- 
ments for two sets of uo n errors. For both control 
procedures, only small differences are noted in the 
flexible-mode responses between the two sets of er- 
rors. As shown in figure 5, the responses null within 
about 10 sec. The relatively small transients in the 
rigid-body responses are primarily caused by the ob- 
server in that perfect estimates of the state variables 
are not possible. (See ref. 4.) In addition, transients 
in the undisturbed modes (initially at 0) are inherent 
to the LQR control procedure. As shown, the rigid- 
body transient responses null in about 60 sec for both 
control procedures. Curves for no-error cases are not 
shown in figure 5 because they would fall within the 
general region of the plotted curves. For larger er- 
rors in uj U) the effects are analyzed in a subsequent 
section of this paper. 

The torque time histories in figure 5 require in- 
stantaneous values because the initial estimates were 
not 0. This condition can be eliminated by incorpo- 
rating small time lags in the control-feedback gain 
matrix. (See ref. 4.) It is of interest to note in 
figure 5(a) that the control actuator requiring the 
maximum absolute value of torque is different for 
the positive- and negative-error cases. This design 
consideration is important in sizing the actuators. 
Figure 5 shows that for LQR control, the total mo- 
mentum for the six actuators is somewhat lower for 
the case with an oj n error of —10 percent. 


The parametric plots in figure 6 show the changes 
in the flexible-mode control requirements from their 
nominal (no-error) values over a range of uj n errors 
of specified magnitudes. The sets for determining 
the random errors used in the frequency-error analy- 
sis are shown in table III. The eight sets represent 
the total number of combinations possible for the 
first three flexible modes. As shown, random errors 
were also incorporated in the three residual modes, 
although it was determined that their effects were rel- 
atively small. The eight sets of random signs shown 
in table III were used to generate constant percent- 
age errors in the six flexible-mode frequencies at spe- 
cific magnitudes of frequency error ranging from 10 to 
50 percent. (It should be noted that in figure 6, data 
for all eight sets of random errors are not shown at 
each specific magnitude of error because some points 
coincide.) The system was stable over the range of 
cu n errors shown, but instability occurred for cu n er- 
rors of ±75 percent. As shown, the effect on total 
momentum is larger for decoupled control for u n er- 
rors of more than 20 percent; the effect on maximum 
torque is larger for LQR control over the entire range 
of t o n error. For both control procedures, it should be 
noted that the 2 -axis control actuator at location 1 
(fig. 1) required the maximum torque in the nominal 
case. In some decoupled cases with an c o n error, the 
2 -axis actuator at location 2 required the maximum 
torque. For the LQR cases represented by the two 
lowest data points across the range of uj u errors, the 
y - axis actuator at location 1 required the maximum 
torque. 

Rigid-body control . The parametric plots in fig- 
ures 7 and 8 show the effects of to n errors on rigid- 
body control requirements. As expected, the changes 
in the control requirements are considerably smaller 
for rigid-body control than for flexible-mode control. 
In figure 7, only the maximum positive and negative 
changes in torque are plotted for the eight cases (ta- 
ble III) at each magnitude of u) n error. From the 
figure, it is evident that effects on maximum torque 
are essentially the same for both control procedures. 
In all cases, the x-axis control actuator at location 2 
required the maximum torque. As for total momen- 
tum requirements, figure 8 shows that u) n errors have 
a negligible effect for decoupled control. 

In reference 4 it was shown that the rigid-body 
modes could be controlled adequately with the use 
of only one three-axis attitude sensor at location 1. 
Hence, an error analysis was made for this condition 
(i.e., no rate sensor). The results are presented 
in table IV for seven sets of random errors in c o n 
ranging from 10 to 20 percent. (No figure is presented 
because variations were relatively small.) The data 


8 



indicate that an unstable control system can result 
with an oj u error greater than 10 percent. For 
example, with an oj n error of 15 percent, the ^ 
rigid-body mode became unstable for LQR control. 
For an uj n error of 20 percent, the same mode was 
unstable along with one of the first three flexible 
modes. For the stable cases, it is seen that the 
maximum torque is unaffected. Also, in two of the 
cases for LQR control with an uj n error of 10 percent, 
the effects on momentum were relatively large as a 
result of poor performance by the £-axis actuators. 
This performance is characterized by the long control 
oscillations that require more than 60 sec to zero 
out, as shown in figure 9. Also shown in figure 9 
is an example of poor response in the ^ rigid-body 
mode and a relatively large transient response for 
the first flexible mode. The </> and 6 rigid-body mode 
responses are not affected. 

Figure 10(a) shows an example of the effect of ob- 
servation spillover (effect of residual modes) on the 
z- axis control actuator at location 1 for decoupled 
control. The other actuators are essentially unaf- 
fected. The rapid control oscillations, although small 
in magnitude, are inherent to the decoupled-control 
procedure, with or without model error, and are ag- 
gravated by the effect of model error. These oscilla- 
tions are eliminated if estimates of the residual modes 
can be included in the control procedure. (See ref. 4.) 
An alternative would be to turn off the actuator, as 
shown in figure 10(b). (This condition was simulated 
on the computer by zeroing out the third row in the 
feedback gain matrix. The other gains in the matrix 
were left unchanged, with the result that the control 
law was not altered.) The system remained stable 
with only minor effects noticed for the first and third 
flexible modes. 

Effect of Model Damping Error on Control 
Performance 

The results in table V show that model damping 
error has little effect on the control requirements, re- 
gardless of whether the error is positive or negative. 
The effect is greater for flexible-mode control; how- 
ever, the larger changes are shown to be negative, a 
result that is not detrimental. (The maximum neg- 
ative error of 50 percent was chosen because some 
natural damping is required for stability.) 

As expected, a negative error has a detrimental 
effect on control actuator performance because of an 
increase in the observation spillover effect caused by 
the reduced damping. The time histories in figure 11 
depict this effect on the performance of the £-axis 
control actuator at location 2. (The other actuators 
are not affected.) Rapid control oscillations for both 
decoupled and LQR flexible-mode control are shown 


in figures 11(a) and 11(b). The LQR oscillations 
in figure 11(b) are more pronounced; however, after 
about 6 sec, their magnitudes are about the same as 
those for decoupled control. (The oscillations can 
be partially alleviated if estimates of the residual 
modes can be included in the control procedure.) 
Figure 11(c) is based on the case for figure 11(b) 
except that the damping error is limited to the first 
three flexible modes. This situation is not likely to 
occur, but the fact that the oscillations are smaller 
than those in figure 11(b) points out the adverse 
effect of negative $ n error on control performance. 
The effect of positive $ n error is not as great, as is seen 
in figure 11(d) for LQR control. Decoupled control 
was not affected. 

Effect of Model Moment-of-Inertia Error on 
Control Performance 

Figure 12 summarizes the effect of model moment- 
of-inertia error on the momentum requirements for 
rigid-body and flexible-mode control. Maximum 
torque requirements are not materially affected. As 
shown, changes in momentum are relatively small 
over the large range of positive and negative errors 
in inertia, especially for flexible-mode control. Note 
that in this case, the effect of LQR control is three 
times larger than for decoupled control. 

Effect of Mode-Slope Error on Control 

Performance 

Consideration of the 15 different nonzero mode 
slopes included in the B and C matrices would yield 
32 768 (i.e., 2 15 ) possible combinations of random 
sets of positive and negative errors for any specified 
magnitude of d>' error. Twenty random sets of er- 
rors (similar to those given in table III) at four val- 
ues of <&' error were selected to determine the max- 
imum changes in control requirements, as shown in 
figure 13. Each random set was selected to be rad- 
ically different from other sets because it was found 
that a change of one or two signs between any two 
sets had a negligible effect on the results; thus each 
set essentially represents 120 or more different sets. 
This value was derived by the following combinato- 
rial formula for the number of r elements out of a set 
of n for a total of 120 combinations. 

c n _ n ( n —!)••• -(re — r+ 1) 
r 1 • 2 • • • -r 

where 

C} 5 = ~ = 15 
and 
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The 20 random sets chosen for the current anal- 
ysis are deemed adequate for establishing the maxi- 
mum effects. The system was stable over the range 
of <&' errors shown. Although errors of ±75 percent 
resulted in a stable system, the control requirements 
were excessive for flexible-mode control. (For exam- 
ple, momentum changes exceeded 300 percent.) 

As expected, the results in figure 13 show only 
small effects of errors on the rigid-body require- 
ments. It is evident that when controlling the flexible 
modes, the maximum torque is more sensitive to 
error for decoupled control than for LQR control. In 
all cases resulting in increased torque for LQR rigid- 
body and flexible-mode control, the control actuator 
requiring the maximum value was the same as that 
for the nominal case. For the decoupled flexible-mode 
control cases resulting in the largest torque increases, 
the z-axis actuator at location 2 required the maxi- 
mum torque. 

Comparison of the results in figure 13 with those 
in figures 6 and 8 shows that the maximum changes 
in momentum are comparable for and uj u errors. 
The larger errors in & are considered to be more re- 
alistic because mode shapes and slopes are the least 
identifiable of the model parameters and are there- 
fore the most important of all the model errors. (See 
ref. 8.) Hence, it would be advisable to concentrate 
on <&' in any initial adaptive control procedure when 
updating the model parameters. 

Effect of Model Error on Maximum 
Transients 

As previously discussed, transients are produced 
in the undisturbed modes when either control pro- 
cedure is employed, but transients are more pro- 
nounced with LQR control. Model errors in uj n and 

affect the maximum values of these transients, as 
summarized in table VI. The last column indicates 
the importance of the particular effect for each con- 
trol procedure. (For example, a small effect repre- 
sents a relatively small change in the transients over 
a large range in model error.) The example in fig- 
ure 14 illustrates the large effect on the flexible mode 
transients due to an oj n error using LQR rigid-body 
control. It is seen in table VI (a) that this effect in- 
creases substantially as the oj n error increases. (Er- 
rors in <; n and in inertia have negligible effects on the 
rigid-body and flexible-mode transients.) 

Figure 15 illustrates a serious effect of uj u and <&' 
errors on the ip rigid-body transient response during 
flexible-mode control for both control procedures. 
Although the maximum transients are only about 
0.2°, it is apparent that a considerable period of time 
is required before the ip responses null. (The other 
two rigid-body responses, 6 and 0, null in 60 sec.) 


For smaller errors in and <&', say 25 percent, the 
ip transient nulls within 60 to 120 sec. 

Effect of Model Error on Time to Null 

Examples of the effect of various model errors 
on the time to null are presented in figures 16 and 
17. Figure 16(a) shows that for large u> n errors 
with decoupled flexible-mode control, the time to 
null is about 60 sec. For the more realistic errors 
in figure 16(c), this time is reduced considerably. 
The results shown in figures 16(e) and 16(f) for large 

errors are similar to those shown in figure 16(a). 
With a few exceptions, model error had no effect on 
time to null for rigid-body control for either control 
procedure. Some LQR cases with oj n errors larger 
than 20 percent required up to 150 sec to null the ip 
response as well as the first flexible-mode transient. 
There is also some effect of model inertia error I 
on time to null, as illustrated in figure 17. (Note 
the expanded time scale.) The curves apply to both 
control procedures. 

Concluding Remarks 

An analysis has been made to determine the ef- 
fects of model error on the control of a large flexi- 
ble space antenna. Control was achieved by employ- 
ing two three-axis control-moment gyros (CMG’s) lo- 
cated on the antenna column. Results have been 
compared for decoupled and linear quadratic reg- 
ulator (LQR) control procedures. Errors were as- 
sumed to exist in the individual model parame- 
ters: modal frequency, modal damping, mode slope, 
and moment of inertia. Their effects on control- 
system performance have been analyzed either for 
(1) nulling initial disturbances in the rigid-body 
modes, or (2) nulling initial disturbances in the first 
three flexible modes. The study includes the effects 
on stability, time to null, and control requirements 
(defined as maximum torque and total momentum), 
as well as on the accuracy of obtaining initial state es- 
timates before the controller is activated. The effects 
on the transients of the undisturbed modes (initially 
at 0) have also been included. 

Errors in the model natural frequencies oj u have 
the largest effect on the accuracy of determining 
initial estimates for the state variables. Model er- 
rors in the mode slopes have the next largest ef- 
fect; errors in the damping and moment-of-inertia 
parameters have a negligible effect. For the wide 
range of model errors investigated, only large errors 
in frequency (> 50 percent) resulted in an unstable 
system. 

Effects of model errors on the control require- 
ments (defined as maximum torque and total 
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momentum) are, in general, comparable for both con- 
trol procedures. The effect of mode-slope error ap- 
pears to be the largest of all model errors. Hence, it 
would be advisable to concentrate on this parameter 
in any initial adaptive control procedure when up- 
dating the model parameters. Errors in the damping 
and moment-of-inertia parameters have only small 
effects on the control requirements. However, nega- 
tive damping error is detrimental to control-actuator 
performance in that rapid oscillations are induced in 
the control responses. 

In some cases for flexible-mode control with fre- 
quency or mode-slope errors, the control actuator 
requiring the maximum absolute value of torque is 
different from the nominal (no-error) case. This dif- 
ference would be an important design consideration 
when sizing the actuators. 


The effects of model errors in frequency and mode 
slope on the maximum transients of the undisturbed 
modes are different for the two control procedures, 
the severity depending upon the type of model error 
and whether rigid-body modes or flexible modes are 
to be controlled. The largest effect is produced 
by frequency errors with rigid-body LQR control. 
Model errors in the damping and moment-of-inertia 
parameters have a negligible effect on the transients. 

For rigid-body control, the model errors have only 
a minor effect on the time to null. Although the 
effect is larger for flexible-mode control, the responses 
nulled within 60 sec in all cases. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
May 27, 1986 
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Appendix 

Tables of Matrices Used in Current Analysis 

This appendix consists of tables of matrices (tables AI to AIV) used in the 
current analysis for the B matrix, F matrix for decoupled control, K matrix 
and Q weightings for LQR control, and the C matrix. 
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TABLE AI. B MATRIX 


B 18x6 


o. 

0. 

0* 

0* 

0. 

0, 

0, 

0, 

0. 

0* 

0. 

0, 

0. 

0. 

0, 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0* 

0, 

0, 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0, 

0, 

0. 

0. 

0, 

0* 

0. 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

•19737917E-07 

0. 

.U003987E-09 

.197379178-07 

0. 

• 1 1003967E-09 

0. 

•19658120E-07 

0, 

0. 

• 196581 20E-07 

0* 

.11003987E-09 

0. 

.25777798E-07 

• 1 1003987E-09 

0. 

•25777798E-07 

.14712325E-05 

0. 

-•1174544 IE— 03 

-.423205428-06 

0. 

-.377761958-06 

0. 

- • 36661 00 3E-03 

0. 

0. 

.31838614E-04 

0, 

.611 43590E-03 

0. 

•89756447E-04 

-.6188551 3 F— 04 

0. 

.20098311E-05 

-•191 7679&E-06 

0. 

•44199044E-05 

.88746034F-07 

0. 

-•10276237E-03 

-# 1 2698521E-03 

0. 

-•11 4521 51E-04 

-• 2 1981 567E-02 

0 . 

-.83552937E-06 

. 59 5 89099E-0 7 

0. 

-•34274029E-05 

• 716896368-07 

0 • 

-.135514 96 E -03 



TABLE AIL F MATRIX USED FOR DECOUPLED CONTROL 




F 6x 12 




Rnw i / -.^67l686lE^05 

Row no. 1 l -. 84090349 E +06 

0« 

. 11921819E +04 

0. 

0 . 

0. 

0, 

.214592758+05 

— .961816728+03 

0, 

-.116065668+04 

Z 1 0. 

-.40617932E+05 

0, 

0. 

0. 

0. 

- • 7 311 22 77E + 06 

0* 

0. 

•22278945E+04 

0. 

• 106861 18E+04 
6 \ • 192 3501 2E + 0 5 

0. 

. 12631355E+04 

0. 

0, 

0. 

0 • 

•22736439E+05 

• 7381 3 68 5E +04 

0* 

-.18777995E+02 

A r - • 459934286+06 
^ -.S278817OE+07 

0. 

•970602748+03 

0. 

0. 

0 . 

0 • 

. 174708498+05 

• 96 181672E +03 

0* 

• 11606566E+04 


-.46807772E+06 

0. 

0. 

0. 

0. 

-.&4253990t+07 

0. 

0. 

-.222789458+04 

0. 

R r . 10<m729E+04 
D '- . 196951 1 2E + 0 5 

0* 

— • 389203 10E +06 

0 • 

0. 

0. 

0. 

-•70056557E+07 

-.73813685E+04 

0. 

• 18777995E+02 


M 

00 



TABLE AIII. K MATRIX and Q WEIGHTINGS USED FOR LQR CONTROL 


Kexl2 


Row no. 1 { 

-.298954E«0t> 

0. 

.54927BE*04 

l643oet^03 

0. 

-.217185E.Q3 

-.496916E+07 

0*. 

. 1 05906 t. 06 

-.379eS2E*03 

0. 

~.26b962t*U4 




2{ 

0. 

-•29J977E* *>6 

0. 

o. 

.2907Q4 L *Q3 

0. 

0, 

•♦493832£*07 

0 . 

0. 

• 3638b9t ^04 

o* 




3{ 

-.643776E.04 

0 * 

- # lbb435£+0b 

* 44b 1 By t ♦ 0 3 

0. 

-.160382E+03 

-.1U494E.06 

0 * 

-.327 1 19E.0 7 

.919151E+04 

0. 

-.3c4594t*U3 




4{ 

-.332541E.06 

0 . 

-.9213ibE>03 

.9iS4iot*00 

0 • 

“.234356E*03 

-.535499E+07 

0* 

13*4 1 1E+05 

,604148E*02 

0. 

.28l73bt*03 




5{ 

0. 

-.337oi2E*Oft 

0. 

0. 

.3l977OE^03 

o. 

0. 

-•543606E* 07 

0. 

0. 

-,342500E>03 

3. 




6{ 

-. 128526E.04 

0 • 

-.205592£*<>6 

.5l*9b5t*03 

0. 

-•221741E.02 

“. 124421E.05 

0. 

-.3B0137E.07 

-.935697E+02 

0. 

-.362612t*0l 




Q 

= diag(0.2 * 

10 12 , 0.2 x 10 12 

, 0.7 x IQ 11 , 

If 1* If 0.2 x 

14 

10 , 0.2 x 10 

14 , 0.92 x io 13 . 

0.41 x 10 8 , 0. 

. 

o 

00 * 
o 

X 


TABLE AIV. C MATRIX FOR THREE- AXIS ATTITUDE-RATE SENSOR AT LOCATION 1 AND 
THREE-AXIS ATTITUDE SENSOR AT LOCATION 4 


C6x18 


Row no. 



22532850F-03 


17983610F-01 
10000000E + 01 
23921470E-07 


37009 500E-04 


0. 

0* 

0. 

0. 

0* 

-.19195410F-02 

0. 

0, 

0. 

0 • 

• 3075645 0F-03 
0. 

aoooooooF+oi 

0. 

0. 

0* 

-.148O60OOE-O4 

0* 


0. 

0. 

.1*790460E-0? 

0, 

0. 

0. 

0. 

0* 

•275R3180F-03 


0. 

- *43030 160E -0 5 


0. 

0* 

0. 

0. 

• LOOOOOOOE ♦01 
-.245M710E-01 
0. 


0* 

• 1 OOOOOOOE + Ol 
-♦1193550QE-06 
0. 

0* 

0. 

0. 

0. 

• 13^647208-05 
-♦1181R08QE-03 
0* 

0. 

0. 

0* 

0* 

-.10959970E-01 

0. 

0. 


0* 

0. 

* * 2 5308020E -03 
0* 

• 10000000F+01 

0* 

0. 

0, 

-•22823800F-04 

0* 

0. 

0 • 

•54079040E-03 

0* 

0. 

0. 

0. 

0. 


0. 

0* 

•43113600E-04 

0. 

0. 

0. 

0. 

•1 OOOOOOOE+O 1 
-.24797850E-02 
-• 5 0625030E-03 

0* 

0. 

0. 

0* 

0. 

.17061950E-03 

0. 

0. 
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TABLE I. ANTENNA CHARACTERISTICS 


'I x = 1.9577 X 10 10 lb-in 2 ; /„ = 1.9657 X lO 10 lb-in 2 ;' 
I z = 1.499 X lO 10 lb-in 2 ; I xz = 0.8357 X 10 8 lb-in 2 ; 
w — 10 020.3 lb 


Mode 

< x ) n , rad/sec 

T n , sec 


Description 

1 

0.7466 

8.42 

0.01 

First torsion 

2 

1.3460 

4.67 

.01 

First bending, zz-plane 

3 

1.7025 

3.69 

.01 

First bending, yz - plane 

4 

3.1813 

1.98 

.01 

Surface torsion 

5 

4.5294 

1.39 

.02 

Second column bending, y^-plane 

6 

5.5905 

1.12 

.02 

Surface-hoop torsion, feed-mast torsion 


TABLE II. SPECIFIED CLOSED-LOOP EIGENVALUES 


Case 

Mode 

rad/sec 

Q 

Decoupled 

0 

0.1 

0.9 



.1 

.9 



,1 

.9 


1 

“.7466 

.59 


2 

“1.346 

.34 


3 

“1.7025 

.24 

LQR 

e 

0.098 

0.904 


<t> 

.112 

.900 



.112 

.901 


i 

.7453 

.502 


2 

1.345 

.511 


3 

1.700 

.502 


“Values used are exact u) n values. 
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TABLE III. RANDOM SETS FOR w„ ERRORS 


Flexible 

mode 

o; n , rad/sec 

Positive and negative errors for random set number — 

i 

2 

3 

4 

5 

6 

7 

8 

1 

0.7466 

+ 


+ 

+ 

— 


— 

— 

2 

1.3460 

+ 


— 

— 

— 

— 

+ 

+ 

3 

1.7025 

+ 

■ 

— 

+ 

— 

+ 

+ 

— 

4 

3.1813 

+ 


+ 

— 

— 

— 

— 

— 

5 

4.5294 

+ 


+ 

— 

— 

— 

— 

+ 

6 

5.5905 

+ 


- 

- 

- 


- 

+ 


TABLE IV. EFFECT OF w n ERRORS ON RIGID-BODY CONTROL REQUIREMENTS 

[ One three-axis attitude sensor ] 


Random error 
in w n , 
percent 

100(T £ —T/y ) /Tjv , 
percent 
(a) 

100(iU € - M n )/M n , 
percent 
(a) 

Comments 

Decoupled 

0 b ) 

LQR 

(c) 

Decoupled 

(b) 

LQR 

(c) 

±10 

0 

0 

1.7 

-1.4 


±10 

0 

0 

-0.6 

10.9 

Long oscillations in 
T Z) \ and T z 2 
for LQR control 

±10 

0 

0 

1.9 

-4.4 


±10 

0 

0 

-0.7 

17.7 

Same as above 

±15 

0 

0 

-0.7 

26.3 

Same as above 

±15 

0 

u 

17.4 

U 

Long oscillations in 
T z , l and T z<2 
for decoupled control 

±20 

u 

u 

U 

U 



a U indicates instability; 0 indicates negligible effect. 

b For decoupled control: T/y = 418 ft-lb; Myy = 10 984 ft-lb-sec. 

c For LQR control: T/y = 301 ft-lb; M/y = 11 720 ft-lb-sec. 
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TABLE V. EFFECT OF j n ERRORS ON CONTROL REQUIREMENTS 


Condition 

Error 
in ? n . 

100(T e - T n )/T n , 
percent 

(a) 

percent 

(«) 


percent 

Decoupled 

LQR 

Decoupled 

LQR 

Rigid-body 

+900 

0 

0.3 

0 

-0.3 

control (one 
three-axis attitude 
sensor) 

-50 

0 

0 

.1 

.1 

Rigid-body 

+900 

-0.7 

0.7 

b 

-0.2 

control (two 
three- axis 
sensors) 

-50 

0 

0 

.2 

0 

Flexible-mode 

+900 

-9.5 

-14.5 

-5.4 

-9.7 

control (two 
three- axis 
sensors) 

-50 

.4 

1.8 

3.7 

4.0 


a Values for T/y and Myy are given in figures 6 to 8 and in table IV; 0 indicates negligible effect. 
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TABLE VI. TYPICAL EFFECTS OF MODEL ERROR ON TRANSIENTS 

(a) Rigid-body control 


Model error 

Error, 

Control 

Maximum transients (absolute) 


in — 

percent 

procedure 

M, in. 

A 2 , in. 

A 3 , in. 

Effect 


±10 

Decoupled 

0.003 

0.014 

0.022 

Small 


±50 


.005 

.040 

.027 



±10 

LQR 

0.21 

0.74 

0.86 

Large 


±50 


.12 

2.17 

2.55 



±10 

Decoupled 

0.002 

0.035 

0.020 

Large 


±50 


.005 

.135 

.175 



±10 

LQR 

0.24 

0.65 

0.75 

Small 


±50 


.21 

.80 

.95 



(b) Flexible-mode control 


Model error 
in— 

Error, 

percent 

Control 

procedure 

Maximum transients (absolute) 

Effect 

0, rad 

<t>, rad 

tp, rad 

0J n 

±10 

±50 

Decoupled 

0.000013 

.000115 

0.000038 

.000137 

0.000405 

.003180 

Large 


±10 

±50 

LQR 

0.000161 

.000194 

0.000233 

.000349 

0.001120 

.001440 

Small 


±10 

±50 

Decoupled 

0.000033 

.000212 

0.000009 

.000146 

0.001040 

.003200 

Large 


±10 

±50 

LQR 

0.000116 

.000078 

0.000270 

.000080 

0.001540 

.003960 

Small 






















ft-lb-sec M vl , ft-lb-sec M„,, ft-lb-sec 


+ 10% errors in w 
-10% errors in w 


1000 

o 

<D 
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Time, sec 

(a) Concluded. 
Figure 5. Continued. 
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Figure 6. Effect of model frequency error on control requirements for decoupled and LQR flexible-mode control. 
The solid and dashed curves were faired through the maximum values obtained. 
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Figure 7. Effect of model frequency error on torque for decoupled and LQR rigid-body control. Only maximum 
positive and negative changes shown; the solid and dashed curves were faired through the maximum values 
obtained. 
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Figure 8. Effect of model frequency error on total momentum for decoupled and LQR rigid-body control. The 
solid and dashed curves were faired through the maximum values obtained. 



ft-lb T„ Q . ft -I 



Time , sec Time , sec 


(a) Torque response. (b) Modal response. 

Figure 9. Example of poor performance for LQR rigid-body control using one three-axis attitude sensor at top 
of antenna column. Error in oj n , ±10 percent; dashed curves indicate no-error case. 
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(a) All six control actuators in operation. 

Figure 10. Typical responses for decoupled rigid-body control using two three-axis sensors. Error in w, 
±10 percent. 
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Figure 13. Effect of model error in <&' on control requirements for decoupled and LQR control. Open symbols 
denote rigid-body control and solid symbols denote flexible-mode control; Tjy and Mpf values are given in 
figures 6, 7, and 8; the solid and dashed curves were faired through the maximum values obtained. 
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